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OF A UNIFORM CANTILEVER BEAM 
By Alexander Mendelscm 


SUMMARY 

Theoretical calculations were made by means of semirigid-flutter 
theory to determine the effect of centrifugal force on the critical 
flutter speed of a uniform cantilever beam. Centrifugal force can 
under certain conditions be detrimental, decreasing the critical 
flutter speed. The higher the ratio of natural torsional freq.uency 
to natural bendisag frequency the greater the decrease in critical 
flutter speed due to centrifugal effect. The conclusion is drawn 
that compressor and turbine blades that are flutter-free at ze 2 ro rpn. 
migh't become unstable at hi^ rotative speeds. The centrifugal force 
effect can therefore not be ignored in any flutter calculations on 
compressor and turbine blades. 


INTRODUCTION 

The vibration problems of compressor and turbine blades have 
become more important as the use of axial-flow compressors and turbines 
has Increased. These problems involve resonant vibrations due to 
mechanical excitation from other parts of the machine or from pulsating 
air flows from previous rows of blades in addition to nonresonant 
aerodynamic excitation. The aerodynamic excitation involves a type of 
vibration that is self-sustained by the continual absorption of energy 
from the air stream. This vibration is called flutter. 

Classical flutter is more specifically defined as a self-sustained 
oscillation due to the coupling of Inertia forces, elastic forces, 
damping forces, and dynamic aerodynamic forces. This type of flutter 
usually occurs on airplane wings at low angles of attack when the 
velocity reaches a certain value called the critical flutter speed. 

It is to be distinguished from a radically different type of flutter 
called stalling flutter, which occurs on profiles at high angles of 
attack, such as propeller blades near the stall point (referaice 1) . 
Stalling flutter can be caused by an aerodynamic hysteresis effect , 
the negative slope of the lift curve, or excitation by a system of 
Karman vortices (references 2 to 4). 
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Classical-flutter theory has been developed by many investigators 
(references 5 to 7) and the results have been found to agree with 
experiment. A direct application of this theory to analysis of com- 
pressor and turbine blades would seem to indicate little lllcelihood 
of such blades fluttering at zero angle of attack. This improbability 
of flutter is due to the stiffness of compressor and turbine blades 
compared to airplane wings. Two Important factors^ however^ must be 
considered; the effect of centrifugal force and the effect of 
cascading. These effects may be either beneficial, that Is, Increase 
the critical flutter speed, or detrimental, that Is, decrease the 
critical flutter speed. 

An approximate method of solutlaa developed at the NACA Cleveland 
laboratory of the effect of centrifugal force on the critical flutter 
speed of a unifona cantilever beam and the results for beams with 
various fundamental freq.uencles are presented herein. The effect of 
cascading is yet to be investigated. 


SYMBOLS 

The following symbols are used in the analysis: 


a 


coordinate of elastic axis measured from midchord In 
units of half chord, positive towards trailing edge 

^ 2 > 

®2 

functions of reduced frequency k and fundamental bending 
frequency (Oj, 


ay, 

a© 

aerodynamic coefficients 

b 


half chord, used as reference unit length 

by. 

by. 

aerodynamic coefficients 

c 


urg2 

e 


pr 

F 


function of k derived by Theodorsen (reference 5) 

fl. 

f2 

functions of ' 
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yo 

P 
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So 


function of k derived by Tbeodorsen. (reference 5) 
moment of inertia about elastic axis per unit span length. 

constant of proportionality 
reduced frequency ~ 

V 

oscil3atory aerodynamic lift force per unit span 
oscillatory aerodynamic moment per unit span 
mass of profile per unit span 

mass of cylinder of air (diameter of cylinder equal to 
profile chord) per unit span 

rotative speed 

location of center cxf gravity of profile measured from 
a in units of half chord 

radius of gyraticm referred to a in units of half chord 

static mass unbalance of profile, mrb 

time 

critical flatter speed 
vertical displacement 
maximum bending amplitude 

phase angle by vhlch torsion vibrational mode lags bending 
mode during flutter 

angle of torsional displacement 

maximum torsional amplitude 

ratio of mass of airfoil per unit span to mass of 

surrounding air cylinder per unit span, “ = — ^- 5 - 

“a 
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P 

(I) 

“b 

time 


density of surrounding air 

frequency of flutter vitration 

fundamental "bending frequency of airfoil 

fundamental "bending frequency of airfoil -corrected for 
centrifugal -force effect 

fundamental torsional frequency of airfoil 

Dots over the symbols represent derivatives with respect to 


ANALYSIS 


Flutter equations. - The analysis of the effect of centrifugal 
force is based on the usual semirigid theory. Reference 8 shows that 
this theory, which assumes two degz^es of freedom, that 1s, the 
uncoupled fundamental bending and torsional modes, gives results in 
close agreement with the exact solution, which assumes an infinite 
number of degrees of freedom. The assumption is made herein that 
these results also apply to the case in which centrifugal forces are 
added. The procedure followed is then to correct the fundam^tal 
bending mode for the centrifugal effect and to proceed in the usual 
manner to obtain the critical speeds and frequencies. The effect of 
centrifugal force on the fundamental torsional mode is assumed to be 
small enough to be neglected. The assumption is also made that the 
oscillatory air forces and momenta acting on the blade are the same 
as those acting on an Isolated two-dimensional airfoil in an uniform 
air stream. Incompressible flow is assumed and structural frictions 
are Ignored. 


The equations for dynamic equilibrium (reference 5) are 
my + S0 + m(Ct3^')^ y - L = 0 


( 1 ) 


10 + Sy + I COt^ 0 - M = 0 
The oscillatory aerodynamic lift L and moment M are given by 
L 


- (ay + lay) y + (a^ + ia@) b0 


( 2 ) 


M = - meCi?b 


(by + Iby) y + (be + Ibe) b© 
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where 



where F and & are functlcaiB of reduced frequencj k given in 
reference 5. 


Let 


then 


itot *1 

y = yo ® I 

e . 9o 


»• p 

y = - 0)'^ y 
0 = - 05^ e 


(4) 


(5) 


By combining eq.uations (l), (2), and (s) and dividing throu^ 
by maCC? 
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y. + lay + 

(^y + ^^y " ®) ^ 


(0‘ 


- 1 


1 


> y + < ag + ’iag - c ^ "be = 0 


bg + ibg + e 


0)2 


- 1 


b^e = 0 


( 6 ) 


Equations (6) have solutions other than y = e o 0 If and only 
if the determinant of the ooefflclents of y and 9 vanishes. By 
setting .this determinant equal to zero, separating real and imaginary 
parts, and rearranging 


and 

\ 

where 

b 


^ 

V“tj " °2 


+ = 0 


( 7 ) 


^ cay + p bg 

C2 = ciy + p bg + ae(by - e) + by(ag - e) - bg ay - ay bg 


a^ = (o - bgKia - ay) + by ag 

bi = cay + bg - c^i ^ 


bg a 


y f®"© “ ^^y 


e) 


. 1 


C]_ = cp 


ej 


For a given profile, a, b, and o are functlonB only of the 
reduced frequency k and the fundamental bending frequency .corrected 
for centrifugal force effects <0|j', which Is, in turn, a functicsti of 
the rotative speed K. For a given rotative speed W, a, b, and o 
are functions only of the reduced frequency k. Equations (7) can 
then be solved graphically for the reduced frequency k and the 
ratio of flutter frequency to fundamental torsional frequency £*)/(o^. 

Bending frequency. - The effect of centrifugal force on funda- 
mental bending frequency C0|, was obtained from reference 9. 

(a>b’)^ = + KW^ 


( 8 ) 
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where cojj is the fundamental tending frequency without rotation 
and £ is a function of the ratio of blade length to root radius and 
of the blade taper. By combining equations (7) and (8), the reduced 
frequency k and flutter frequency W can be obtained as a functlcai 
of the rotative speed N. The critical flutter speed v can then be 
obtained from the relation 


V 


(Ob 

k 


( 9 ) 


HESOTJTS AND DISCUSSION 

Equaticais (7) were solved for the reduced frequency k and the 
frequency ratio as functions of the ratio of torsional frequency 

to bending frequency corrected for centrifugal-force effect ' .* 

the results are shown in figure 1. A typical profile with the following 
constants was chosen for the calculatlcaiB : 

a = - 0.40 

r = 0.10 
rg = 0.50 

p = 100 

Fran figure 1 and equation (8) it is possible to obtain the 
reduced frequency k and the frequency ratio co/io^ as a function of 
the rotative speed N, once a value for the fundamental biding 
frequency ccjj is chosen. For this purpose, equation (8) can be 
conveniently rewritten in the form 



For a given profile, the frequency ratio CO^/Co^’ will therefore 
decrease as the rotative speed N increases. F!r*om figure 1, for 
values of frequency ratio less than approximately 0.90, the 

value of reduced frequency k will decrease and the value of (O /6)^ 
will decrease as the rotative speed N increases (or as 
decreases) . At values of frequency ratio ' greater than 

approximately 0.90, both reduced frequency k and frequency 
ratio increase as the rotative speed N Increases (or as 

decreases). But for values of close to 0.90, the 



8 


MACA BM No. E8B05 


reduced frequency k Increases at a faster rate than 
critical flutter speed v can he expressed as 

V = — 

k 

The critical flutter speed t will therefore increase with increasing 
rotational speed N for values of frequency ratio ' less 

than 0.90 and will decrease with increasing rotative speed N for 
values of co^/co^j' greater than 0.90. The existence of a minimum 
value for the critical flutter speed v is therefore indicated. 

The foregoing illations can perhaps he more clearly seen if 
equation (lO) is used with figure 1 to plot directly the reduced 
frequency k and the frequency ratio w/O-t aa a function of the 
rotative speed N. This variation is shown in figures 2 and 3 for a 
fundamental bending frequency of 100 radians per second and 

frequency-ratio ' values ranging from 1 to. 5. The constant Z 

/ 2 rt \2 

was assumed to he • This value correspcHids to a non- 

• tapered blade with a ratio of root radius to blade length of zero. 

As shown in figures 2 and 3, the frequency ratio Increases 

as the rotative speed N Increases but the reduced frequency k 
at first increases and then decreases. If the rate of increase of 
reduced frequency k in any part of the region where it is 
increasing is greater than the rate of Increase of the frequency 
ratio CO/4o^, then from equation (ll) a minimum critical flutter 
speed V can be expected. 

In order to verify this conclusion, the critical flutter 
speed V was calculated assuming a blade chord of 2 inches. The 
results are shown in figure 4 for values of fundamental bending 
frequency 0>b of 100, 500, 1000, and 2000 radians per second, 
respectively. The critical flutter speed v decreases to a minimum 
value as the rotative speed N is increased. This minimum value 
depends on the bending frequency cOjj, on the ratio of torsional to 
bending frequencies and on the rotative speed N. The 

rotative speed N at which the critical flutter speed v is a 
minimum Increases as the ratio of torsional to bending frequency 

increases and as the natural bending frequency increases; 

also, the higher this ratio, the greater the decrease in critical 
flutter speed v. This decrease in critical flutter speed v as 
the rotative speed W is increased is due to the centrifugal force 
causing only the bending frequency co-j^ to increase. Although 
increasing both the bending frequency cOj, and torsional 
frequency will Increase the critical flutter speed v, increasing 


The 
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the bending freq.uency 03|j alone may decrease the critical flutter 
speed (reference 10) . This relation can he seen if the curves of 
figure 1 are expressed as 


^ \^'J 

^ 

Then hy eq_uaticai (ll) 

■ ^ '■© 


If the torsional freq,uency 
increased so that the ratio 


remains constant^ 



co^ and bending freq.uency ' are both 
of torsional to bending freq.uency (O-fc/o-jj' 

will remain constant. The value of 


critical flutter speed v, however. Increases because of the inci*ease 
in torsional frequency 0)^. A change in the value of bending 
frequency OJ-jj' alcne, however, may either raise or lower the value 



speed V. 


and produce 


a corresponding change in critical flutter 


At high frequency ratios and at high values of the 

natural bending frequency w-jj, the critical flutter speeds v 
obtained are above the lower compressibility limit for which the cal- 
culations are valid and the results are therefore meaningless. This 
limit has been arbitrarily chosen as 800 feet per second and all 
curves above that speed are shown as deished lines. 


For natural bending frequencies COf, hi^er them 1000 radians per 
second, classical flutter cannot occur in the subsonic range except at 
low frequency ratios OOfc/cOjj and possibly at high rotative speeds N 
as shown by figure 4(c) and 4(d). Because current compressor and 
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turbine blades usually have bending frequencies greater than 
1000 radians per second, little danger of- such blades fluttering 
classically would exist if the torsional frequency (0^ were several 
times the bending frequency However, at high rotative speeds N, 

the critical flutter speed v can be greatly reduced and blades that 
would be safe at zero rpm can become unstable. As an example, fig- 
ure 4(c) shows that at a frequency ratio of 2 the flutter speed v 

is reduced from a value that is out of the subsonic range to about 
615 feet per second as the rotative speed is increased frcxQ 0 to 
10,000 rpm. 


CONCLUSIONS 

From the calculations made to determine the effect of centrifugal 
force on the flutter of a uniform cantilever beam, the conclusion 1s 
drawn that under certain conditions centrifugal force can be detri- 
mental, decreasing the critical flutter speed. A compressor or tur- 
bine blade that would be stable at zero rpm might become unstable at 
high rotative speeds. The flutter characteristics of blades /oust 
therefore be investigated at the operating engine speeds. The 
validity of the semirigid flutter theory with the centrifugal force 
effect included must still be verified. The effect of cascading must 
also be investigated. 


Flight Propulsion Eesearch Laboratory, 

National Advisory Comaittee for Aeronautics, 
Cleveland, Ohio. 
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Figure 1, - Variation of ratio of flutter frequency to 
torsional frequency and of reduced frequency with 
ratio of torsional frequency to corrected bending 
frequency. 
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to bending frequencies • Fundamental bending frequency, 100 
radians per second. 
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0 2,000 4,000 6,000 8,000 10,000 

Rotative speed, N, rpm 

(a) Fundamental bending frequency, 100 
radians per second. 

Figure 4. - variation of critical flutter speed with, 
rotative speed for various ratios of torsional to 
bending frequencies. 
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0 2^300 4.000 6,000 8,000 10,000 12,000 14,000 16,000 


RotatlTO speed, K, rpa 

(b) Fundamental banding frequency, 600 radians per second. 



Figure 4. - Continued, variation of critical flutter speed with rotative speed for various 
ratios of torsional to bending frequencies. 


critical flutter 



Rotative speed, K, rpn 


(c) Pvindamental bending frequency, 1000 radians 
per second. 


Figure 4. - Continued, Variation of critical flutter speed with, 
rotative speed for various ratios of torsional to bending 
frequencies. 








Critical flutter speed, v, ft/seo 


18 


NACA RM No. E8B05 


1800 


1600 


1400 


1800 


1000 


BOO 


600 


400 


200 


0 



2,000 4,000 6,000 8,000 10,000 12,000 
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(d) F^mdanientsl bending frequency, 2000 radians per second. 

Figure 4. - Concluded. Variation of orltlcal flutter speed ultli rotative speed for various 
ratios of torsional to bending frequencies. 







